function [] = fig_sec10_4_congestionControlDual()

close all;

hd = load ('data/congestionControlDual_hd.txt');
pie = load ('data/congestionControlDual_pie.txt');
ye = load ('data/congestionControlDual_ye.txt');
objFunc = load ('data/congestionControlDual_objFunc.txt');
jom_hd = load ('data/congestionControlDual_jom_hd.txt');
jom_pie = load ('data/congestionControlDual_jom_pie.txt');
jom_objFunc = load ('data/congestionControlDual_jom_objFunc.txt');

multiplot ('Time', 'Network Utility' , objFunc , jom_objFunc , [-11 0] , '-' , '\pi_e' , pie , jom_pie , [0 0.05] , '-' , 'h_d' , hd , jom_hd , [0 30] , '-');

% multiplot ('Time', 'Network utility' , objFunc , jom_objFunc , [-900 -500] , '-' , 'p_e' , pe , jom_pe , [0 0.2] , '-' , 'u_e' , ue , jom_ue , [0 0.1] , '-');
% multiplot ('Time', 'Network utility' , objFunc , jom_objFunc , [-100 100] , '-' , 'p_e' , pe , jom_pe , [0 1] , '-' , 'u_e' , ue , jom_ue , [0 1] , '-');

jom_objFunc
grad_objFunc = objFunc (end)

max(hd(end,2:end))
max(jom_hd)

